Exo-hydrogenated Single Wall Carbon Nanotubes 
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An extensive first-principles study of fully exo-hydrogenated zigzag (n,0) and armchair (n,n) single 
wall carbon nanotubes (CnHn), polyhedral molecules including cubane, dodecahedrane, and CeoHgo 
points to crucial differences in the electronic and atomic structures relevant to hydrogen storage 
and device applications. CnH„'s are estimated to be stable up to the radius of a (8,8) nanotube, 
with binding energies proportional to f/R. Attaching a single hydrogen to any nanotube is always 
exothermic. Hydrogenation of zigzag nanotubes is found to be more likely than armchair nanotubes 
with similar radius. Our findings may have important implications for selective functionalization 
and finding a way of separating similar radius nanotubes from each other. 



PACS numbers; 61.48. +c,61.46.-fw,61.50.Ah,71.15.-m,71.20.Tx 



Carbon nanotubes 0] exhibit very unusual structural 
and electronic properties, suggesting a wide variety of 
technological applications [§|j3|, including the storage of 
hydrogen where the large effective surface area promises 
a large absorption capacity P-^10|]. Unfortunately, the 
studies to date report conflicting results. While some 
labs [^,|| report hydrogen storage densities up to 10 wt% 
other labs report only 0.4 wt% on the same system. 
Theories based on physisorption have failed to predict 
such high uptake ||ll[]. To the best of our knowledge, 
studies of hydrogen chemisorption in nanotubes are very 
limited [ p"2[p^ ] and are clearly needed to have a better 
understanding of hydrogen and nanotube system. 




Fig. 1 Three different polyhedra of carbon and hydrogen; 
(a) Cubane (Oh), (b) Dodecahedrane (Ih), and (c) a side and 
top view of a single wall exo-hydrogenated carbon nanotube. 

Hydrogen-carbon interactions have been studied ex- 
tensively both theoretically and experimentally for many 
interesting polyhedral molecules, such as cubane (CgHg) 
[p^-p^, dodecahedrane (C20H20) 0, and various iso- 
mers of CgoHn [^,|9| (see Fig. 1). Despite its very 
strained 90° CCC-bond angle cubane has been synthe- 



sized successfully ||lj] (Fig. la). Similarly, dodecahedrane 
and various isomers of CgoHn (up to n=32) have been 



also synthesized |18|. These novel polyhedral molecules 
which represent the zero dimensional case, exhibit many 
interesting properties. However, due to the one dimen- 
sional nature and the curvature of carbon nanotubes, 
the hydrogen-carbon interactions in these systems may 
be quite different than those in polyhedral molecules. 
Therefore, it is important to know if it is also possible to 
hydrogenate carbon nanotubes in a similar way and if so 
what their structural and electronic properties would be. 
This paper addresses this important issue by perform- 
ing extensive first-principles calculations and shows that 
the chemisorption of hydrogen is dependent on the radius 
and chirality of the nanotubes. Theoretical predictions 
from first-principles studies played an important role in 
guiding experimental studies in the past [ p^j20| and we 
expect that many findings reported here may have im- 
portant implications in this interesting system as well. 

In order to obtain a reasonably complete understand- 
ing, we studied a very large number of systems including 
zigzag (n,0) (n=7,8,9, 10 and 12) and armchair (m,m) 
(m=4,5,6,8, and 10) nanotubes |gj as well as cubane, do- 
decahedrane, CgoHeo and finally hydrogenated graphene 
sheet (i.e. infinite limit of tube radius). The first prin- 
ciples total energy and electronic structure calculations 
were carried out using the pseudopotential plane wave 
method ||2^. The results have been obtained within the 
generalized gradient approximation (GGA) This 
method has already been applied to many carbon sys- 
tems, including fuUerenes and cubane with remarkable 
success jl^,^ . We used plane waves with an energy cut- 
off of 500 eV. With this cutoff and using ultra soft pseu- 
dopotentials for carbon j2j], the total energy converges 
within 0.5 meV/atom. Interactions between molecules or 
nanotubes in periodic cells are avoided by using large su- 
percells. The supercell parameters are chosen such that 
the closest H-H distance is 6 A. For molecular calcu- 
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lations the Brillouin zone integration is carried out at 
the F-point. For nanotubes, we used k-point spacing 
of dk « 0.02 A"-'^, generating 5 and 10 special k-points 
along the tube-axis for zigzag and armchair nanotubes, 
respectively |2^ . All carbon and hydrogen positions were 
relaxed without assuming any symmetry. For nanotube 
calculations, the c-axis of the supercell (corresponding to 
the tube axis) is also optimized. 
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Fig. 2 Energy curves as a CH-bond is rotated towards 
the indicated arrows for armchair (top) and zigzag nanotubes 
(bottom) starting from 9 = a. The minimum energy is found 
when the CH-bond is tilted toward the shaded hexagons. Zero 
of energy was taken to be arbitrary. 

In principle, there are an infinite number of isomers 
depending on the locations of hydrogen atoms (i.e. endo 
if they are inside the tube and exo if they are out- 
side) as well as amount of hydrogen coverage. Endo- 
hydrogenation, alternating endo-exo hydrogenation, and 
various half coverage cases are being studied and the re- 
sults will be published elsewhere |2^. Here, we consider 
the case of full coverage where all carbon atoms in a nan- 
otube are hybridized with hydrogen atoms from outside 
of the tube as shown in Fig. Ic. We refer to this isomer 
as a fully exo-hydrogenated carbon nanotube. 

First, the equilibrium orientations of the CH bonds 
were determined starting with all the CH-bonds radially 



outward (Fig Ic). Using this configuration, we studied 
a single CH-bond orientational dependence of the poten- 
tial energy surface. Figure 2 shows the calculated energy 
curves as a single CH-bond is rotated along two high sym- 
metry directions for both zigzag and armchair nanotubes. 
For armchair (n,n) nanotubes, the optimum orientation 
is obtained when the CH-bond is tilted about the tube 
axis (i.e c-axis). Hence in the fully optimized structure, 
CH-bonds tilt in opposite directions around the c-axis al- 
ternatively. For the zigzag nanotubes, the optimum ori- 
entation is obtained when the CH-bond is tilted towards 
the c-axis. Therefore, the lowest energy configuration 
for zigzag tubes has CH-bonds tilted towards plus and 
minus c-axis alternatively. Having located the CH-bond 
orientations in this way, we next let all the carbon and 
hydrogen atoms along with c-axis vary to obtain the final 
optimum structures. 

Table 1 summarizes the parameters obtained for fully 
optimized structures. Upon hybridization of carbons 
with hydrogens, the C-C bond length {dec) increases 
from « 1.4 A to « 1.55 A. The latter is typical for sp'^ 
CC-bonds. The increase in dec results in an increase 
in the tube radius (Rhc) by about 13 - 16 % for arm- 
chair nanotubes and by about 15 - 17% for zigzag nan- 
otubes. Interestingly, these values are almost twice of 
those found for the polyhedral molecules. Moreover, the 
value of dec increases slightly (by about 0.03 A) with in- 
creasing tube radius. The CH-bond length {dcH « 1.09 
A) is also found to have weak dependence on the tube 
radius. Using projection techniques we estimated the to- 
tal charge transfer from hydrogen to carbon to be around 
0.26 electrons for nanotubes and 0.3 electrons for poly- 
hedral molecules. 

The most important difference between zigzag and 
armchair nanotubes is found in the local CCH angles 
{c(cch)- Even though one of these angles is about the 
same for both types of nanotubes, the second angle in 
zigzag nanotubes is always larger than those in arm- 
chair nanotubes. This implies that the CCH-bond an- 
gles are more frustrated in armchair nanotubes than in 
zigzag nanotubes and therefore deviate more from the 
ideal tetrahedral sp'^ bond angle of 109.5°. This observa- 
tion suggests that hydrogenated armchair nanotubes will 
have higher energy and therefore they are less stable than 
zigzag nanotubes. Unlike CCH-bond angles, CCC-bond 
angles have weak radius dependence and are about the 
same for both types of nanotubes. 
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Material 


Formula 


Rhc ( Rc) (A) 


dcH (A.) 


■ dec (A) 


accH (deg.) 


accc (deg.) 


(4,4) 
(5,5) 
(6,6) 

(8,8) 
(10,10) 


C20H20 
C24H24 
C32H32 
C40H40 


3.103 (2.734) 
3.885 (3.394) 
4.698 (4.061) 
6.228 (5.400) 
7.780 (6.755) 


1.090 
1.087 
1.084 
1.079 
1.077 


1.541, 1.567 
1.549, 1.575 
1.557, 1.594 
1.567, 1.594 
1.574, 1.600 


96.70, 98.60 
94.82, 97.15 
93.35, 96.30 
92.16, 94.85 
91.40, 94.00 


112.77, 120.69 
113.18, 121.50 
113.30, 122.00 
114.62, 121.95 
115.40, 121.76 


(7,0) 
(8,0) 
(9,0) 
(10,0) 
(12,0) 


C32H32 
CieHig 

C40H40 

C48H4g 


3.180 (2.765) 
3.641 (3.146) 
4.111 (3.557) 
4.571 (3.912) 
5.467 (4.695) 


1.092 
1.090 
1.089 
1.087 
1.084 


1.549, 1.553 
1.553, 1.557 
1.553, 1.566 

1.556, 1.572 

1.557, 1.576 


96.40, 102.25 
95.22, 101.60 
94.32, 101.14 
93.60, 100.54 
92.66, 99.340 


113.95, 125.90 
114.12, 127.00 
114.27, 127.58 
114.48, 127.85 
115.11, 127.67 


Cubane 
Dodecahedrane 
FuUerene 


CsHg 
C20H20 


1.345 (1.267) 
2.157 ( 2.0) 
3.827 ( 3.510) 


1.087 
1.090 
1.088 


1.553 
1.539 
1.536 1.550 


125.26 
110.9 
101.1, 101.9 


90.0 
108.0 
108.1, 120.0 



Table 1 Various parameters of the fully optimized structures of exo-hydrogenated armchair and zigzag carbon nanotubes 
and other polyhedral molecules. For graphene (i.e. Rhc ^ 00) dcH and dec are 1.066 A and 1.622 A, respectively. 



The stability and energetics of CH-bond formation are 
derived from the average binding energy per atom for 
exo-hydrogenated nanotubes defined as 

Eb = {Ec^H„ -Ec- nEH)/n. (1) 

Here Ec„h„j Ec, and Eh are the total energies of the 
fully optimized exo-hydrogenated nanotube, nanotube 
alone, and hydrogen atom, respectively. According to 
this definition, a stable system will have a negative bind- 
ing energy. Fig. 3a shows the radius dependence of Eb 
for nanotubes and polyhedral molecules (see inset). Two 
interesting observations are apparent. First, as shown by 
solid and dotted lines, the binding energies can be very 
well described by a one parameter fit; 

EB=Eo-C{n,m)/RHc, (2) 

where Eq is the limit Rhc ~^ 00 (i.e. graphene) and 
calculated to be -1.727 eV. The fit results for C{n, m) are 
given in Fig. 3a for zigzag and armchair nanotubes. The 
inset to Fig. 3a shows that while for cubane falls on 
the same curve as nanotubes, dodecahedrane and CeoHgo 
have lower energies than nanotubes due to the their more 
spherical shape. 

The second interesting observation in Fig. 3a is that 
the binding energies of zigzag nanotubes are always lower 
than those in armchair nanotubes with similar radius 
by about 30 meV/atom. As discussed above, this is a 
natural result of the fact that the CCH-bond angles in 
zigzag nanotubes are closer to the optimum tetrahedral 
sp^ bonding than those in armchair nanotubes. We ex- 
pect this observation is also valid for hybridization of 
nanotubes with other elements, such as CI and F and this 
may have important implications for separating similar 
radius nanotubes from each other by selective chemical 
functionalization. 
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Fig. 3 (a) Binding energies Eb of (n,n) (square) and (n,0) 
(circle) nanotubes as a function oiYinc- The solid and dashed 
lines are one parameter fits to Eb — Eq — C{n,m)/ Rhc as 
discussed in the text. Inset shows the binding energies of 
cubane, dodecahedrane, and CeoHeo. (b) Full circles indi- 
cate the energy (AEb) to break a single CH-bond to form a 
CnHn-i zigzag nanotubes as depicted in the top inset. Full 
squares indicate energy gain AEq by attaching a single H 
atom to a nanotube to form a CnH as depicted in the bottom 
inset. The solid and dashed lines are two parameter fits as 
discussed in the text, indicating 1/Rhc behavior. 
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Even though C„H„ nanotubes are found to be sta- 
ble with respect to a pure carbon nanotube (C„) and 
n X H atoms for all values of the radius, it is of interest 
to see if they are also stable against breaking a single 
CH-bond. We, therefore, calculated energies of fully op- 
timized hydrogenated nanotubes after breaking one of 
the CH-bonds and putting the H atom at the center of 
supercell as shown in the top inset to Fig. 3b. Calcu- 
lated values of the energy differences AEb, for zigzag 
nanotubes were fitted to AEb = Eq + A/Rhc where Eq 
and A are 2.506 eV and -15.671 eVA, respectively. We 
note that for radius around Rhc ~ 6.25 A, the AEb be- 
comes negative, suggesting instability Hence, (12,0) 
and (8,8) nanotubes are at the limit for stable, fully exo- 
hydrogenated nanotubes. We are currently studying this 
problem for half-coverage case as well. 
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Fig. 4 (a) Electronic band structure of a (9,0) exo- 
hydrogenated carbon nanotube (left panel) and the corre- 
sponding density of states (DOS) (right panel) . (b) Band gap 
as a function of tube radius Khc- Inset shows the full scale 
plot to include the band gaps of the polyhedral molecules. 

The energy, AEq, gained by attaching a single H atom 
to a carbon nanotube is calculated by performing struc- 
ture optimization of a C„H-nanotube as depicted in the 
bottom inset to Fig 3b. It is seen that AEq can be also 
well described by AEq — E'q + A' /Rhc (dashed line) 



where E'g and A' are -1.161 eV and -4.952 1 eVA, re- 
spectively. Unlike AEb, there is no change in the sign 
of AEg, suggesting that for any radius of carbon nan- 
otube hybridization of a single carbon atom is always 
stable. However the energy gain from two such processes 
is around 5-6 eV which is slightly less than the dissosia- 
tion energy of H2, 6.65 eV. Hence C„ nanotube plus H2 
system is stable against forming a C„H2 hydrogenated 
nanotube. Therefore, in order to realize the CH-bonding 
discussed here, one first has to break H2 molecules into 
hydrogen atoms, probably by using a metal catalyst or 
electrochemical techniques. 

Hydrogenation of nanotubes is also important in the 
modification of the electronic structure for device ap- 
plications. Figure 4 shows the band structure and the 
corresponding density of states (DOS) for a (9,0) exo- 
hydrogenated nanotube, which is typical to other nan- 
otubes that we studied. Using projected DOS, we find 
that the bottom of the conduction bands are mainly de- 
rived from hydrogen while the top of the valence bands 
are mainly carbon-origin. In contrast to pure nanotubes 
which are metal or semiconductors depending on their 
structure, the C„H„ nanotubes are found to be direct 
band insulators with a gap of 1.5-2 eV at the F-point. 
This value is about one-third of those for the molecu- 
lar polyhedrals, indicating less stability of hydrogenated 
nanotubes than molecules (Fig. 4b). The band gaps de- 
crease with increasing tube radius but unlike binding en- 
ergies there is no apparent I/Hhc type behavior. Inter- 
estingly, the band gaps of armchair nanotubes are higher 
in energy by about 0.2 eV than those in zigzag nanotubes. 
This is surprising because the band gap is usually higher 
for more stable saturated hydrocarbons. 

The observed band gap opening via hydrogenation of 
nanotubes can be used for band gap engineering for de- 
vice applications such as metal-insulator heterojunctions. 
For example, various quantum structures can easily be 
realized on an individual carbon nanotube, and the prop- 
erties of these structures can be controlled by partial hy- 
drogenation of carbon nanotubes. If the different regions 
of a SWNT are covered with hydrogen atoms, the band 
gap and hence the electronic structure will vary along 
the axis of the tube. This way various quantum struc- 
tures of the desired size and electronic character can be 
formed. In this respect, present scheme is quite similar to 
our previous constructions of nanotube heterostructures 
or quantum dots, where periodic applied transverse com- 
pressive stress is used for band gap opening p^ ]. 

In summary, we have presented first-principles calcula- 
tions of the structural and electronic properties of various 
nanotubes which arc fully protonated by sp^ hybridiza- 
tion of carbons. Wc find that C„H„ nanotubes are stable 
for tube radius Kmc smaller than 6.25 A, roughly corre- 
sponding to a (8,8) nanotube. Hybridization of a single 
carbon atom is found to be always exothermic regardless 
of tube radius. Weak but stable CH-bonding in nan- 
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otubes may be an important consideration for possible 
hydrogen storage applications. We also found that hy- 
bridization of zigzag nanotubes is more likely than arm- 
chair nanotubcs with the same radius, suggesting a pos- 
sible selective chemical functionalization of nanotubes. 
The fact that other carbon clusters such as cubane, do- 
dccahcdranc, and C60H32 have been synthesized success- 
fully, suggest that it may possible in the near future to 
hydrogenate carbon nanotubes, yielding new structures 
with novel properties. 

ACKNOWLEDGEMENTS. This work is partially 
supported by the National Science Foundation under 
Grant No. INT97-31014 and TUBITAK under Grant 
No. TBAG-1668(197 T 116). 



[1] S. lijima, Nature 354, 56 (1991); S. lijima, T. Ichihashi, 

and Y. Ando, Nature 356, 776 (1992). 
[2] N. Hamada, S. Sawada, and A. Oshiyama, Phys. Rev. 

Lett. 68, 1579 (1992). 
[3] M. S. Dresselhaus, G. Dresselhaus, and P. C. Eklund, 

Science of Fullerenes and Carbon Nanotubes; Academic 

Press, San Diego (1996). 
[4] A. C. Dillon, K. M. Jones, T. A. Bekkadahl, C. H. Kiang, 

D. S. Bethune, and M. J. Heben, Nature 386, 377 (1997). 
[5] C. Liu, Y. Y. Fan, M. Liu, H. T. Cong, H. M. Cheng, 

and M. S. Dresselhaus, Science 286, 1127 (1999). 
[6] C. M. Brown, T. Yildirim, D. A. Neumann, M. H. Heben, 

T. Gannett, A. C. Dillon, J. L. AUeman, and J. E. Fis- 
cher, Chem. Phys. Lett. 329, 311 (2000). 
[7] C. C. Ahn, Y. Ye, B. V. Ratnakumar, C. Witham, R. C. 

Bowman, and B. Fultz, Appl. Phys. Lett. 73. .378 (1998). 
[8] P. Chen, X. Wu, J. Lin, K. L. Tan, Science 285, 91 

(1999). 

[9] Qinyu Wang, S. R. Challa, D. S. ShoU, and J. K. Johnson, 

Phys. Rev. Lett. 82, 956 (1999) 
[10] R. C. Gordillo, J. Boronat, and J. CasuUeras, Phys. Rev. 

Lett. 85, 2348 (2000). 
[11] V. V. Simonyan, P. Diep, and J. K. Johnson, J. Chem. 

Phys. Ill, 9778 (1999). 
[12] K. Tada, S. Furuya, and K. Watanabe, Phys. Rev. B 63, 

155405 (2001). 

[13] Y. Ma, Y. Xia, M. Zhao, R. Wang, and L. Mei, Phys. 

Rev. B 63, 115422 (2001). 
[14] P. E. Eaton and T. W. Cole, Jr., J. Am. Chem. Soc. 86, 

962 (1964). 

[15] T. Yildirim, P. M. Gehring, D. A. Neumann, P. E. Eaton, 
and T. Emrick, Phys. Rev. Lett. 78, 4938 (1997). 

[16] T. Yildirim, S. Ciraci, and A. Buldum, Phys. Rev. B 62, 
7625 (2000). 

[17] B.S. Hudson, D. A. Braden, S. F. Parker, and H. 

Prinzbach, Angew Chem. Int. Edit. 39, 514 (2000). 
[18] M. I. Attalla, A. M. Vassallo, B. N. Tattam, J. V. Hanna, 

J. Phys. Chem. 97, 6329 (1993). 
[19] T. Guo and G. E. Scuseria, Chem. Phys. Lett. 191, 527 



(1992). 

[20] T. Yildirim, O. Giilseren, Q. Kilig, and S. Ciraci, Phys. 
Rev. B 62, 12648 (2000). 

[21] Structure of a nanotube is described by two integers 
(n,m) as explained in M. S. Dresselhaus, G. Dresselhaus, 
and R. Saito, Phys. Rev. B 45, 6234 (1992). 

[22] M. C. Payne, M. P. Teter, D. C. Allen, T. A. Arias, and J. 
D. Joannopoulos, Rev. Mod. Phys. 64, 1045 (1992). The 
computer code CASTEP is distributed and maintained 
by Molecular Similations Inc. Identification of commer- 
tial products does not imply recommendation or endorse- 
ment by the National Institute of Standards and Tech- 
nology. 

[23] J. P. Perdew and Y. Wang, Phys. Rev. B 46, 6671 (1992). 
[24] D. Vanderbilt, Phys. Rev. B 41, 7892 (1990). 
[25] H. J. Monkhorst and J. D. Pack, Phys. Rev. B 13, 5188 
(1976). 

[26] O. Giilseren, T. Yildirim, and S. Ciraci, in prep. (2001). 

[27] There is an energy barrier for breaking a single CH-bond 
process from the fully hydrogenated carbon nanotube, 
AEb is the energy difference between the initial and final 
states. 

[28] Q. Kihg, S. Ciraci, O. Giilseren, and T. Yildirim, Phys. 
Rev. B 62, R16345 (2000). 



5 



